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ABSTRACT 

Based on recent models of relativistic jet formation by thermal energy depo- 
sition around black hole-torus systems, the relation between the on- and off-axis 
appearance of short, hard gamma-ray bursts (GRBs) is discussed in terms of 
energetics, duration, average Lorentz factor, and probability of observation, as- 
suming that the central engines are remnants of binary neutron star (NS+NS) 
or neutron star-black hole (NS+BH) mergers. As a consequence of the interac- 
tion with the torus matter at the jet basis and the subsequent expansion of the 
jets into an extremely low-density environment, the collimated ultrarelativistic 
outflows possess flat core profiles with only little variation of radially-averaged 
properties, and are bounded by very steep lateral edges. Owing to the rapid 
decrease of the isotropic-equivalent energy near the jet edges, the probability 
of observing the lateral, lower Lorentz factor wings is significantly reduced and 
most short GRBs should be seen with on-axis-like properties. Taking into account 
cosmological and viewing angle effects, theoretical predictions are made for the 
short-GRB distributions with redshift z, fluence, and isotropic-equivalent energy. 
The observational data for short bursts with determined redshifts are found to 
be compatible with the predictions only if either the intrinsic GRB rate density 
drops rapidly at z > 1, or a large number of events at z > 1 are missed, implying 
that the subenergetic GRB 050509b was an extremely rare low-fluence event with 
detectable photon flux only because of its proximity and shortness. It appears 
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unlikely that GRB 050509b can be explained as an off- axis event. The detection 
of short GRBs with small Lorentz factors is statistically disfavored, suggesting 
a possible reason for the absence of soft short bursts in the duration-hardness 
diagram. 

Subject headings: gamma-ray bursts — hydrodynamics 

1. Introduction 

The origin of short GRBs has long been the subject of theoretical speculation. A long- 
standing prediction is that short GRBs might originate from NS+NS/BH mergers (e.g., 
Blinnikov et al. 1984; Paczyhski 1986; Eichler et al. 1989; Meszaros 2002, and references 
therein). In the widely favored scenario a NS+NS/BH system gives rise to a BH girded by 
a dense torus of gas from the disrupted NS companion. The BH accretes mass at rates up 
to many solar masses per second, releasing huge amounts of gravitational binding energy 
mostly in neutrinos, gravitational waves, and through magnet ohydro dynamic processes. 

Hydrodynamic simulations (Aloy, Janka & Miiller 2005; AJM05) show that thermal en- 
ergy deposition around post-merger BH-torus systems, e.g., by neutrino- antineutrino annihi- 
lation, can drive collimated, ultrarelativistic outflows with the high Lorentz factors, internal 
variability and internal shocks that are deemed necessary to explain GRBs with the fireball 
model (e.g., Piran 2005). The collimation of observed short GRBs is rather uncertain and a 
matter of ongoing studies. Evidence for collimation was claimed in case of GRB 050709 (Fox 
et al. 2005), for which a jet semi-opening angle of 9- ]Ct ~ 14° was inferred. Although more 
observations are needed to confirm the collimation of short GRBs, this first indication seems 
to agree with the hydrodynamic simulations of AJM05. If the majority of short bursts has 
opening angles of that size, it would mean that only about 1% of all bursts point to Earth. 

Only in the past year have a few short GRBs been localised. Their properties are 
recapped in Table 1. Two of these short GRBs are hosted by elliptical galaxies, and also the 
host galaxy of GRB 051221 appears to have a relatively evolved population of stars, despite 
its higher star-formation rate. This adds confidence to the idea that the progenitors of short 
GRBs are old systems. In addition, GRB 050509b is in the outskirts of its potential host, 
in agreement with expectations that coalescing compact binaries can travel large distances 
away from their birth sites during their long gravitational-wave driven inspiral (Tutukov & 
Yungelson 1994; Bloom, Sigurdsson, & Pols 1999). Moreover, a possible SN association has 
been ruled out for GRB 050509b (Hjorth et al. 2005; Castro-Tirado et al. 2005) and for 
GRB 050709 (Covino et al. 2005). 
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Three out of five short GRBs were detected at similar redshifts (z ~ 0.2), but their 
intrinsic isotropic-equivalent energy release in gamma rays, E^^, spans a factor of ~ 100. 
The observation of intrinsically less energetic short GRBs at higher redshifts may be selected 
against. 

In this paper we use the jet models of AJM05 to discuss theoretical predictions for the 
observable properties of short GRBs as a function of viewing angle (Sect. 2). In Sect. 3 we 
investigate the consequences of these models for the probability of observing short bursts 
with different isotropic equivalent energies and Lorentz factors, including the selection effects 
due to cosmological redshift. A summary and conclusions follow in Sect. 4. 

2. Jets from post-merger black holes 

By means of relativistic hydrodynamics simulations AJM05 studied the formation of 
ultrarelativistic outflows from BH-torus systems. The BH and torus masses and the depo- 
sition of thermal energy around the BH were chosen as expected from NS+NS/BH merger 
models. For different assumptions about the environment density and for varied geometry 
and time-dependence of the energy deposition rate, AJM05 followed the acceleration, colli- 
mation, and propagation of the ultrarelativistic outflows for an evolution time of 0.5 seconds. 
If the deposition rate of thermal energy per solid angle around the BH was sufficiently large, 
ultrarelativistic jets were launched along the rotation axis 1 . Although these simulations are 
still far from taking into account all potentially relevant physics and did not self-consistently 
track the viscosity-driven torus evolution and its neutrino emission, they nevertheless provide 
useful insight into the conditions and properties of mass outflows from post-merger BH-torus 
systems. 

The fact that an ultrarelativistic jet is launched from a "naked" BH-torus system has 
important consequences for its properties. While in the case of collapsars as sources of 
long GRBs the ultrarelativistic outflow originates from BH-torus systems at the center of 
a massive star, the polar jets in NS+NS/BH mergers do not have to plough through many 
solar masses of overlying stellar matter. Since the acceleration is not damped by swept-up 
matter, the jets very quickly reach Lorentz factors of a few. The acceleration is driven by 
the enormous radiation pressure of the pair-photon fireball produced by the energy release 
around the BH. Collimation of the baryon-poor jets is provided by the much denser torus 



1 A lower limit to the energy deposition rate needed to drive a relativistic outflow is set by the fact that 
the ram pressure of polar mass infall to the BH must be overcome (see AJM05 for details; model B03 in that 
paper did not reach the energy deposition threshold and is therefore not considered here). 
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walls which gird the evacuated polar regions of the BH. As the jets propagate into the 
extremely low-density environment, they continue to accelerate, reaching maximum Lorentz 
factors of a few 100 within the simulated evolution time. 

As a consequence of the interaction with the torus matter at the jet basis and the 
subsequent free expansion, the collimated ultrarelativistic outflows possess flat core profiles 
with only little variation of radially-averaged specific properties. These cores are bounded 
at their lateral edges by very steep gradients, which are not smoothed by the prolonged 
entrainment of baryons as in case of long GRB jets. The rapid decrease of the isotropic- 
equivalent energy as a function of 9 implies, among other effects, that the probability of 
observing the lateral, lower Lorentz factor wings must be expected to be significantly reduced. 

Figure 1 provides an overview of the hydrodynamics results for the sample of "type- 
B" models discussed by AJM05. In these simulations the BH-torus system was assumed 
to be surrounded by a very rarified medium. The models differ in the adopted rate of 
thermal energy deposition per solid angle and in the prescribed time dependence. The 
parameters were chosen such that the neutrino-antineutrino annihilation as calculated from 
NS+NS/BH merger simulations, and in particular for the post-merger accretion of a BH 
(Ruffert & Janka 1999; Janka et al. 1999; Janka & Ruffert 2002; Setiawan, Ruffert & Janka 
2004), was qualitatively and quantitatively reproduced. Because of large variations of the 
torus masses expected from compact object mergers 2 , AJM05 explored a variety of energy 
deposition properties, approximately covering the range of predictions from the above merger 
simulations (cf. Table 1 in AJM05). 

The corresponding differences in the assumed energy deposition around the BH-torus 
systems account for the variation of the on- axis values of the isotropic- equivalent kinetic 
plus internal energies, E^iY^), of jet matter with estimated terminal Lorentz factors > 
100 (Fig. 1, left), is estimated by assuming that 50% of the total specific energy will 
ultimately be converted to kinetic energy at large radii 3 . The models span a range of ~ 100 in 
.Eiso, from some 10 49 erg to several 10 51 erg. Depending on the efficiency of energy conversion 
to 7-ray emission, the measurable 7-ray energy (E 1;iso ) may be roughly a factor of 10 lower 
than E^T^) (e.g., Guetta, Spada & Waxman 2001, and references therein). All models 
show a flat core and steep lateral wings in the radial average of V (Fig. 1, right), but models 



2 The mass of the remnant torus varies with the masses and spins of the binary components and with 
their mass ratio; it is also sensitive to still uncertain properties of dense NS matter and depends on general 
relativity effects that arc included only in a subset of the current NS+NS/BH merger simulations. 

3 Note that an exact calculation of is beyond the scope of the hydrodynamic models of AJM05, because 
the simulations do not include physics that plays a role during the later stages of the jet acceleration. 
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B07 and B08 have a somewhat more shallow decline of E iso because a gradual decrease of 
the energy deposition rate was assumed after an initial, burst-like phase. This caused a 
continuous decrease of the jet opening angle, thus softening the wings of the jet profile. It is 
currently not clear whether the jet-driving energy release of the BH-torus system follows such 
a burst/slow-decay behavior or whether it is powerful and roughly constant over a longer 
period of time before it ceases more abruptly, as assumed in the other models (B1-B6). 

Figure 2 provides some observationally relevant quantities deduced from the hydrody- 
namics results. The profiles of ^^(Too) given in the left panel there represent the isotropic- 
equivalent values of energy which is potentially radiated from the outflow in different di- 
rections. These profiles are calculated from the jet energies of Fig. 1, taking into account 
radiation contributions coming from regions outside of the line of sight. In order to compute 
these contributions, let us consider an amount of energy dE' that is radiated into the solid 
angle dfl' at angle 9' relative to the direction of motion (primed quantities are measured in 
the local comoving frame). The transformation of the energy per solid angle between the 
comoving frame and the laboratory frame (e.g., a frame at rest with respect to the central 
BH) is (Rybicki & Lightman 1985) 



where (3 is the velocity of the comoving frame as measured in the laboratory frame. Com- 
paring the values of ^ at two different angles, 9 = 9 1 and 9 = 0, one finds that the energy 
emitted per unit of solid angle around an angle 9\ can be expressed in terms of the energy 
radiated along the direction of propagation (9 = 0) as 



if the radiation field is isotropic in the comoving frame. The emission from the jet in a 
certain observer direction 9\ can now be obtained as the sum of the contributions from all 
matter moving in different directions with velocities /3o at angles 9 relative to the jet axis, 
yielding 



Since the profiles of E iso (9) displayed in Fig. 1 are related to the energy radiated per unit solid 
angle into different directions 9 = 9 , the integration given in Eq. (1) can be directly applied 
to these profiles in order to obtain the isotropic equivalent energy distributions corrected for 
off-line-of-sight contributions. The results are plotted in Fig. 2. From the latter figure it 
is obvious that energy release with -Ei SO (roo) > 10 48 erg is confined to semi-opening angles 
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between about 10° and 20°, corresponding to a collimation ("beaming") factor of the two 
polar jets of fn = 1 — cos#j et between 1.5% and 6%. 

The relative observability of a short burst from a model of our set at different polar 
angles 9 is shown in Fig. 2, middle panel. We plotted there the quantity 

_ sin 9 E iso (9) 
j:\dcos0' E iso (9') ' 

where as a rough estimate for the detectability of a burst we assumed fdet(9) oc ^^(Too) 
(using for ^^(Too) the data from the left panel in Fig. 2). The probability distribution 
confirms the steep lateral edges visible in the other quantities for models B01-B06 and 
the slightly softer wings in models B07 and B08, which correlate with a sharp drop of the 
average Lorentz factor. The probability P{9) peaks at angles between 4° and ~11°, which 
is a somewhat smaller range than that associated with £ , i so (r oo ) > 10 48 erg. Morever, P{9) 
nicely demonstrates the widening of the lateral visibility of a GRB with increasing energy 
deposition per solid angle around the BH (cf. the sequence of models B02, B04, B06, B01, 
and B05). 

The hydrodynamic jet models of AJM05 track the formation of ultrarelativistic GRB- 
viable outflow. It turns out that the "shells" ejected by the central engine, i.e. the BH-torus 
system, accelerate much faster in the leading part of the outflow than the shells in its lagging 
part. The rear shells need therefore a longer time to reach velocities v « c. This different 
acceleration at early and late times of the relativistic wind ejection leads to a stretching of the 
overall radial length of the outflow, A, relative to the on-time t ce of the central engine 4 times 
the speed of light c, A > ct cc (see Fig. 3). This stretching has the important consequence 
that the overall observable duration of the GRB (in the source frame), T — tA — A/c, may 
be a factor of ten or more longer than the activity time of the central energy source (Fig. 3), 
even when the GRB is produced by internal shocks. This is only seemingly in conflict with 
the canonical fireball picture as discussed, e.g., by Piran (2005); Kobayashi, Piran & Sari 
(1997); Sari & Piran (1997); Nakar & Piran (2002), who draw the conclusion that "internal 
shocks continue as long as the source is active, thus the overall GRB duration T reflects 
the time that the inner engine is active". It should in fact be noted that GRB light curve 
calculations in these shell collision models are based on the assumption that the shells move 
with constant Lorentz factors and with a velocity very close to the speed of light, which means 
that the expansion of the shells is considered only after the preceding acceleration away from 



4 We define the on-time or activity time of the central engine as the period of time during which the 
energy release of the BH-torus system is sufficiently powerful to drive an ultrarelativistic outflow as required 
for GRBs. 
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the central engine. Therefore the "shell emission time from the inner engine" (Piran 2005) or 
"ejection time" (Kobayashi, Piran & Sari 1997) or "activity time of the inner engine" (Nakar 
& Piran 2002) referred to in these works corresponds to the overall length the fireball has 
attained in the saturated stage after the initial acceleration. This time can only be identified 
with the real activity time or on-time of the central engine, whose energy output launches 
and drives the initial shell ejection, if the shells accelerate to nearly the speed of light within 
a negligibly small period of time. Only in this case the claim of equal emission and source 
activity times applies and can be traced back to the fact that a photon radiated from a 
fireball shell is observed almost simultaneously with a (hypothetical) photon launched at the 
central engine together with the radiating shell (Nakar & Piran 2002). The hydrodynamic 
jet simulations, however, show that the acceleration time, in particular for the shells ejected 
later, is not negligible (Fig. 3). 

Taking the terminal length A(9) of the ultrarelativistic outflow in direction 9 as a very 
crude measure of the overall observed duration T(8) of a burst 5 , we have plotted this duration 
relative to the on-axis value in the right-hand panel of Fig. 2. The off-axis variations of the 
possible GRB durations are sizable. Bursts which are seen more than 5°-10° off-axis have 
only 30%-60% of the on-axis duration. The outflow stretching can easily be the dominant 
effect determining the duration of short GRBs, but it is not possible to assess its relevance 
for long bursts without hydrodynamical modeling. This is due to the more complex and 
largely different dynamics of the collimated relativistic outflow caused by the presence of the 
massive star in the collapsar model. 

Predictions of the absolute duration of the GRB emission based on jet simulations where 
the period of energy release was just a parameter are hard to make. Self-consistent torus 
evolution models, including the energy production mechanism and the feedback from jet 
formation, will be needed to provide reliable numbers for the duration of the source activity 
and the differential acceleration of the ejecta shells at early and late times. Moreover, the jet 
structure at 0.5 s after the jet creation, at which time AJM05 had to stop their calculations, 
may not be representative for the situation at the time when the GRB is produced. 



5 We estimate A(#) as the radial extension of those parts of the outflow which reach terminal Lorentz 
factors > 100, where of a mass element in the outflow is again calculated by assuming that 50% of 
its total energy will finally be converted to kinetic energy. 
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3. Predicted short GRB distributions 

We can employ our set of jet models to attempt to predict the observable short-GRB 
distributions versus redshift z, fluence, isotropic equivalent energy output in ^y-rays, Ey t i SCn 
and Lorentz factor T of the flow that produces the bulk of the measured 7-radiation. We 
discuss these aspects in view of the observed bursts of Table 1. Our methodology is different 
from that of other groups, who used observational data of GRB redshift, luminosity, and peak 
flux distributions to derive constraints on the intrinsic properties (event rates or lifetimes, 
jet collimation, luminosities) as functions of redshift (e.g., Ando 2004; Guetta & Piran 2005; 
Nakar, Gal- Yam & Fox 2005). We instead refer to the sample of jet models from AJM05 
to define the intrinsic distribution of short GRB energies as a function of viewing angle for 
the individual events. This approach suffers from the huge uncertainties of the theoretical 
models (see below) and naturally can have only a tentative character. 

In order to perform the analysis, we assume that our models have equal probability 
and are representative of the intrinsic variability of the short GRB source population at all 
redshifts. This is a crucial assumption and it can certainly be questioned, but currently it 
can hardly be replaced by more realistic alternatives, because the uncertainties concerning 
the link between NS+NS/BH binary parameters and GRB properties are still large. 

In order to test the sensitivity of our analysis to a variation of the comoving short GRB 
rate density as function of redshift, Rqrb(z), we consider three cases, namely an intrinsic 
GRB rate which (i) is constant with comoving volume, (ii) follows the star formation rate 
SF2— sfr of Eq. (4) in Guetta & Piran (2005), or (iii) varies according to the binary neutron 
star merger rate, case SF2+delay of Eq. (6) in Guetta & Piran (2005). 

Based on these prescriptions, we produce the expected observable distributions of GRB 
properties by Monte Carlo sampling, randomly drawing GRB energies and Lorentz factors 
as functions of viewing angle 6 from our set of GRB-jet models, using the results shown in 
Fig. 2 and Fig. 1, respectively. We estimate the isotropic-equivalent GRB energy output, 
E 7> i so (0), by reducing the isotropic equivalent kinetic plus internal energy of the outflow, 
Ei SO (9) of Fig. 2, by a factor of 10. This is supposed to account for the limited efficiency of 
energy conversion to 7-rays and the fact that only a fraction of the 7-emission occurs in the 
energy band of a measurement 6 . Our set of models can be considered to define a co-moving 
space distribution function E lt i so , F), which is assumed to be normalized to unity and 
whose integral over F yields another normalized distribution function, <&(9, E lt i so ) = J dF $. 



6 A careful inclusion of a redshift-dependent k-correction appears inappropriate, because we do not apply 
detailed theoretical arguments to estimate the conversion efficiency of jet energy to 7-radiation, and model- 
dependent spectral properties of the 7-ray emission are disregarded as well. 
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The random angle between the jet axis and the line of sight is denoted by 6. Formally, 
the expected redshift distribution of measured bursts can be computed from the intrinsic 
distribution as 

+1 

N( Zl ,z 2 ) = jdz^^f^jd^n J dE*{0,E), (3) 

when N(zi, z 2 ) represents the observed rate of events in the redshift interval z\ < z < z 2 (in 
the following we will only consider normalized distributions and the absolute value of the 
rate will be of no concern). In Eq. (3), \i = cos9, dV/dz = 4nD 2 L (z)c[H (l + z) 2 (Q M (l + z) 3 + 
^(I+^+^a) 1 ^ 2 ] -1 is the comoving volume element, and (l+z) _1 accounts for cosmological 
time dilation. The cosmological parameters used in our study are H = 72kms~ 1 Mpc~ 1 , 
tt A = 0.72, tt M = 0.28, and Q K = 1 - - fi A . In E( l- ( 3 ) as wel1 as in Ec l s - ( 4 H 6 ) 
below, an instrument specific detection probability that depends on the photon flux and 
thus on the burst luminosity, spectrum, and redshift, is ignored (we have no information 
on this quantity in terms of the fluence). i?™™ o (0, / m in, z) is the isotropic equivalent energy 
corresponding to the minimum fluence / m i n that can be measured by the detector, and 
E™?^(9) corresponds to the maximum energy release of the models of our sample in a certain 
direction 9. The fluence at the instrument is computed from the intrinsic energy output 
E lt i so (d) (in the reference frame of the source) and from the luminosity distance Dl(z) as 
f(9,z) = (l + z)E J - lso (0)/[^Dl(z)]. 

The predicted distribution of Lorentz factors of the observed bursts can be written as 

+ 1 T2 ^max(-E'y,iso(^))/min) 

N(ri,r 2 ) = Jd^ Jdv J d£ 7iiso $(#,£ 7)iso ,r) J dz^- RG ^ ] , (4) 

-1 Ti 

where z rnax {E 1 ^ so {9), / m i n ) is the maximum redshift at which GRBs with isotropic equivalent 
energy E lt i so (9) produce a fluence above the lower detection bound / m i n . Similarly, the 
fluence distribution is given by 

+ 1 OO £ 7 ,iso(/2 0)) 

N(f\,f 2 ) = Jd^nj dz^^lM J dE<S>(9,E) (5) 

- 1 £ 7 , iso (/i(z)) 

with f 2 > fi > /min, when E 1;iso (fi(z)) and E 1 ^ so {f 2 {z)) are the isotropic equivalent energies 
of 7-radiation that account for fluences f\ and f 2 in the frequency window of the detector 
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for a GRB at redshift z. The distribution versus E^ :iso can be represented by 

+ 1 "^7,iso ^ ; max(-E'7,iso(^))/min) 

iV« iso ,< iso ) = |d /U 2 7 r J d£ 7 , iso *(0,£ 7 , iso ) J ^^TVT- (6) 



E 1 . 

"y,iso 



We perform our investigation with two different values of the lower fluence cutoff for 
detection (i.e., the detection threshold is assumed to be a step function). For one simulation 
we use / min = 10~ 8 erg cm~ 2 , which corresponds to the fluence of the very weak GRB 050509b 
and thus can be considered as a lower bound of the sensitivity of Swift. In a second simulation 
we use a cutoff value of f m - m = 1.6 x 10~ 7 ergcm~ 2 , which we derive from the limiting flux, 
0min = /min/[(l + ^T^], of I photon cm~ 2 s~ 1 adopted for the BATSE instrument by Guetta 
& Piran (2005), making the assumption that {l J r z)T- l e 1 = lOOkeVs is a representative value 
for the product of photon detection time interval and energy (TJ is the intrinsic duration of 
the burst). This higher fluence threshold seems to be compatible with the sample of (bright) 
BATSE bursts listed in Table 1 of Ghirlanda, Ghisellini & Celotti (2004). We note that Nakar 
et al. (2005) employ a detection threshold of 1 photon cm" 2 s~ 1 also for Swift. A higher value 
for /min may further be motivated by the fact that the photon flux decreases with increasing 
redshift more rapidly than the fluence by an additional factor (I + z)~ l . Therefore a low 
fluence threshold for detection may overestimate the number of high- ,2 events observed. 

To obtain normalized distributions, N(xi,x 2 )/N, we do not evaluate the integrals in 
Eqs. (3)-(6) directly, but perform Monte Carlo sampling of a large number of GRB events at 
different redshifts, with random orientations and with properties drawn randomly from the 
model data plotted in Figs. 1 and 2, thus constructing the discrete distribution functions $ 
and $ which describe the GRB properties according to our set of jet simulations. The sam- 
pled events with / > / min are then collected into bins to build up the normalized probability 
distributions versus redshift, fluence, E^^, and V shown in Fig. 4. We point out that the 
Lorentz factors are based on those of Fig. 1 and therefore reflect the situation only 0.5 s after 
the ejection of relativistically expanding matter. Because of subsequent acceleration and 
conversion of internal to kinetic energy (which we are unable to trace in the hydrodynamic 
simulations), the terminal values must be expected to be larger by a factor of 2-3. 

Our predictions for the observable distributions vs. redshift, fluence, energy, and Lorentz 
factor for the two chosen values of the fluence cutoff are displayed in Fig. 4. The correspond- 
ing locations of four of the five bursts listed in Table 1 are also indicated. The very bright 
short GRB 051221 is outside of the energy scale of the third panel, hence we omitted it 
from Fig. 4. This burst is too energetic to be accounted for by the range of 7-burst energies 
considered in our set of jet models, in which we assumed that 10% of the ejecta energy can 
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be converted to radiation in the frequency band of the measurement. But if this fraction 
were significantly higher — the observations of GRB 051221 indeed suggest a total efficiency 
of 60-70% (Soderberg et al. 2006) — then the output of 7-energy of this burst would be 
within the reach of our models for near-axis observation (cf. Fig. 1). Our jet models in- 
voked energy deposition rates as obtained for the annihilation of neutrino-antineutrino pairs 
in NS+NS/BH merger and post-merger accretion simulations (for details and references to 
original work, see AJM05). We therefore conclude that GRB 051221, if it is indeed lo- 
cated at a redshift of 0.546, does not make a strong case for different underlying energy 
extraction mechanisms (e.g., magnetohydrodynamic processes) and/or different progenitors 
as speculated by Soderberg et al. (2006). 

It is obvious that the choice of the lower fluence cutoff has a big impact on the redshift 
distribution (left panels). For a detection threshold of 10~ 8 ergcm~ 2 and our employed co- 
moving GRB rate densities, the majority of bursts should be seen at redshifts z > 1, while 
with the higher cutoff we predict that bursts with z > 0.75 would not be detected. Clearly, 
the observations of the recent bursts with redshifts between 0.16 and 0.72 are more compat- 
ible with the latter case. At the request of one referee, we performed a Kolmogorov-Smirnov 
(KS) test (not including the very energetic GRB 051221), which gives significance levels of 
5%, 8%, and 20% for the consistency of the observations with the theoretical distributions for 
the three tested GRB rate densities SF2— sfr, SF2+delay, and "constant", respectively (the 
KS probabilities in case of the lower cutoff value are 0.1%, 0.4% and 2%). Better agreement 
is therefore found when the simulation predicts a larger number of nearby events. Highest 
preference is attributed to the constant GRB rate density, in which the relative fraction of 
visible bursts at z < 0.5 is largest. This is in qualitative agreement with the recent analysis 
by Nakar et al. (2005; see Fig. 5 there), who found that the well- localized bursts imply a 
high local rate of short-hard GRBs events, and that long lifetimes of the progenitor systems 
are favored. Given the small number of events, the results clearly will change significantly 
with every new observation (see Bloom & Prochaska 2006). 

The fluence distributions show the expected increase towards faint events, and the distri- 
bution of -E 7i iso reveals a strong bias towards events with higher isotropic-equivalent energies, 
corresponding to near-axis observation. Such events are much more likely to yield a fluence 
above the threshold values at large redshifts and are also preferred because of the steepness of 
the jet wings, which reduces the probability of seeing events from those wings. GRB 050509b 
is an exceptionally weak event, which is well separated from the other observed bursts, and 
which the theoretical distributions do not account for. It thus distorts the KS measures 
for the fluence and energy distributions, which otherwise signal reasonably good compati- 
bility between calculations and observations, in particular for the higher fluence threshold. 
GRB 050509b is most probably not an off-axis observation, but requires either an event 
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with lower energy output than for any of the jet models in our sample, or an efficiency of 
the energy conversion to 7-rays lower than assumed here. The latter is disfavored by the 
observations of Bloom et al. (2005). These authors also discard GRB 050509b as an off-axis 
event on the basis of the very early decay of the afterglow light curve. 

The probability distribution of the Lorentz factor reflects the discreteness of our sam- 
ple of models, because the redshift has no influence on T and because the jet wings are 
extremely steep and the off-axis visibility much reduced (see Fig. 1). Therefore, the high- 
energy model B05 with a relatively wide jet clearly sticks out in the T-distribution. A 
continuous set of models would enhance and broaden this peak, and the "underwood" of 
events extending to low values of T would even be more reduced if terminal Lorentz factors 

instead of T were plotted. 

We hypothesize that this clear bias towards observing high-r events may be connected 
with the lack of short-soft GRBs in the duration-hardness diagram. Short bursts are on 
average harder than long bursts, and short-soft bursts are rare or missing, because short 
GRBs develop jets with higher Lorentz factors and very steep edges. The underlying reason 
behind this effect is the fact that jets from post-merger BH-torus systems do not have to 
make their way out of a massive star through layers of dense stellar material, which can 
damp the outflow acceleration by mass entrainment. In the commonly used internal shock 
scenario for the GRB emission the peak energy of the (synchrotron) spectrum, E p , obeys the 
relation E p oc T b with b ~ —2 (e.g., Daigne & Mochkovitch 2003; Ramirez-Ruiz & Lloyd- 
Ronning 2002; Zhang & Meszaros 2002). In view of our jet models this would imply that 
short bursts might have lower E p but nevertheless be harder than long ones because of a 
steeper increase of the vF v spectrum below the peak. The results of Ghirlanda et al. (2004) 
seem to confirm this. Their spectral analysis for a sample of short bright GRBs detected 
by BATSE in comparison with the spectral properties of long bright BATSE bursts indeed 
reveals a significantly harder (i.e., less negative) spectral index a rather than a larger peak 
energy E p for short GRBs. This leads to an increased hardness ratio if E p is above the 
higher energy band of the spectral hardness measure, or overlaps with it, which is the case 
for almost all of the short GRBs listed in Table 1 of Ghirlanda et al. (2004). We point out 
that the difference between the peak spectral energies of short and long bursts might even be 
more pronounced than found by Ghirlanda et al. (2004) if long bursts occurred at typically 
higher redshifts and the intrinsic values of these energies were compared. 

However, the discussion by Zhang & Meszaros (2002) also leaves possibilities how in the 
internal shock model higher Lorentz factors of short GRBs may be compatible with higher 
peak spectral energies. In case of Poynting-flux dominated outflow, Zhang & Meszaros 
(2002) expect a direct proportionality, E p oc T (see their Eq. 21 and case (b) in their 
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Fig. 3), and for a kinetic-energy dominated fireball they find E p oc L 1 / 2 /?^ (see their Eq. f 7) 
when R int T 2 cSt is the radius where internal shock collisions produce 7-ray emission (see, 
e.g., Piran 2005). Although short GRBs have higher T's compared with long GRBs, their 
variability timescale St might be generically much smaller. This would result in a smaller 
shock dissipation radius and therefore in a higher E p , because magnetic fields may be stronger 
at a smaller radius. 



4. Discussion and Conclusions 

Using a set of hydrodynamic simulations of ultrarelativistic jet formation by thermal 
energy release around BH-torus systems as remnants of NS+NS/BH mergers (AJM05), we 
have investigated the off-axis visibility of short GRBs produced by such outflows. The jets are 
characterized by narrow cores with high isotropic-equivalent energies and very high Lorentz 
factors, which are laterally bounded by steep wings. These properties are a consequence of 
the specific conditions in which the jets are considered to be launched, i.e., their accelera- 
tion away from the torus-girded BH, where the energy is released, into an environment of 
extremely low density. The observability of short GRBs is therefore strongly favored within 
a cone of semi-opening angle between about 10° and 15° around the jet axis. 

We argued that the observable duration of short GRBs both in the external shock 
model and in the internal shock collision model for GRB production can be significantly 
longer than the activity time of the central engine whose energy release powers the GRB- 
viable ultrarelativistic outflow. The reason for this claim is the fact that our jet formation 
simulations showed that during the acceleration phase the fireball experiences significant 
radial stretching, because hydrodynamic effects cause the acceleration to proceed differently 
in the leading and lagging parts of the outflow. Therefore the overall radial length A of the 
relativistic ejecta can become a factor of ten or more larger than ct ce , when t ce is the on-time 
of the central engine (Fig. 3). This seems to be in conflict with the canonical internal shock 
scenario, according to which the observed GRB is produced by shell collisions such that 
the observed light curve reflects the temporal activity and overall duration of the energy 
release by the "inner engine" (Piran 2005). However, the standard picture applies only if the 
fireball shells move immediately after their ejection with velocities that are very similar (but 
not necessarily identical) and very close to the speed of light. This can be traced back to 
the fact that a photon radiated from a fireball shell is observed almost simultaneously with 
a (hypothetical) photon emitted from the central engine together with the radiating shell 
(Nakar & Piran 2002). If the acceleration is differential, and in particular if the acceleration of 
shells ejected by the central engine at later times proceeds more slowly, the fireball evolution 
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after creation at the engine is more complex (Fig. 3). Its emission pattern then does not 
directly replicate the activity timescales and the on-time of the central source. 

In order to investigate the consequences of the reduced off-axis visibility of short GRBs 
produced by our jets, we made an attempt to predict observable short-hard GRB distribu- 
tions with redshift, fluence, isotropic-equivalent energy, and Lorentz factor. For this purpose 
we made use of a sample of seven ultrarelativistic jet models (AJM05), which roughly covers 
the range of energy release expected from neutrino-antineutrino annihilation in NS+NS/BH 
binary mergers. Employing these models for defining the intrinsic GRB properties, we per- 
formed a Monte Carlo sampling of the expected observed event distributions. Our analysis 
took into account the off-axis variation of the jet properties as provided by the hydrodynamic 
simulations (Figs. 1 and 2) and referred to standard prescriptions for the intrinsic event rate 
density as a function of redshift (constant rate as well as the cases SF2— sfr and SF2+delay 
from Guetta & Piran 2005). We found that the resulting distributions predict far too many 
bursts at z > 1 and therefore show no good agreement with the recent observations when 
a fluence threshold of 10 -8 ergcm~ 2 as suggested by GRB 050509b was adopted (the KS 
probability for the observed and predicted data being drawn from the same distribution is 
less than 1%). Improved consistency was achieved when we performed our analysis with a 
higher fluence threshold of 1.6 x 10~ 7 ergcm -2 . In this case the bursts at z > 0.8 escape 
from observation and the KS significance increases to ~ 20% for the constant comoving GRB 
rate density. This would imply that GRB 050509b was a very special and extremely rare 
case with an exceptionally low fluence, which due to its shortness and proximity still had a 
sufficiently large photon flux of about 1 photon s" 1 cm' 2 and was therefore above the Swift 
detection threshold. Comparison of the results for different comoving GRB rate densities 
revealed that the consistency between theory and observations increases when the rate yields 
a relatively larger number of events at redshifts z < 0.5 rather than an increase for higher 
values of z. Our results agree qualitatively with those of Nakar et al. (2005), who found that 
the detection of the current sample of \ow-z short GRBs is best compatible with a very low 
comoving rate density for z > 1 and a peak at z < 0.5. 

The extremely weak GRB 050509b is clearly separated in energy from the other ob- 
served events and cannot be accounted for by off-axis observation of any of the modelled jet 
outflows. It seems to require an intrinsically less energetic event, because an efficiency of 
energy conversion to 7-rays significantly lower than what we used is disfavored by observa- 
tions (Bloom et al. 2005). On the other hand, the very energetic GRB 051221 exceeds the 
isotropic-equivalent GRB energies predicted by our jet models, making the assumption that 
10% of the ejecta energy can be converted to radiation in the detector frequency band. If, 
however, this fraction is as high as suggested by the measurements (60-70% for the total 
efficiency in case of GRB 051221; Soderberg et al. 2006) and GRB 051221 happened indeed 
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at a redshift of 0.546 and not farther away, then this strong burst is also within the reach of 
the most energetic one of our models if observed near-axis. Therefore we do not think that 
GRB 051221 gives strong support to speculations that different progenitors and/or different 
energy extraction mechanisms than neutrino-antineutrino annihilation are needed to explain 
the observed large spread in short burst energies and the high energy of GRB 051221 in 
particular (Soderberg et al. 2006). 

Because of the steep wings of the jet profiles in both isotropic-equivalent energy and 
Lorentz factor V, our analysis reveals a clear bias towards the detection of short GRBs with 
high T values. Since higher Lorentz factors and steep jet edges may be characteristic features 
of GRB jets that originate from post-merger BH-torus systems, which makes them different 
from collapsar jets, we propose that they are the reason why short GRBs are typically harder 
than long ones. In the internal shock scenario for GRB emission, this might imply that the 
peak energy of the synchrotron spectrum is lower (E p cx R^ n l oc r~ 2 5t -1 ). Short bursts could 
therefore be harder not because of a higher E v but because of a steeper increase of their vF v 
spectra below the peak energy. This would require the spectral maximum to lie above or 
inside the high-energy band for the hardness ratio. Such properties were indeed concluded 
from a spectral analysis of a sample of short bright BATSE GRBs in comparison with long 
burst spectra (Ghirlanda et al. 2004). Alternatively, a higher T could still allow for a larger 
E p if the variability timescale 5t of short GRBs were generically much smaller than that of 
long GRBs. 

We point out that our exploration can only be tentative. Both the theoretical models 
and observational data still involve large uncertainties, which have an impact on our analysis. 
Not only did we take a very simplistic approach in estimating the 7-ray production of our 
jet models in the energy band of the detector. We just assumed that a fraction of 10% 
of the ejecta energy is converted to 7-emission at the relevant frequencies. We also did 
not account in detail for the instrument specific detection probability that depends on the 
photon flux and thus on the burst luminosity, spectrum, and redshift, but simply assumed 
that all short GRBs above a certain fluence threshold are detected, independent of their 
peak luminosity and spectral properties. Moreover, the intrinsic short GRB rate per unit 
of comoving time and comoving volume enters the analysis sensitively. Nakar et al. (2005) 
found that the well-localized short GRBs, which are all at z < 1, require a very large number 
of nearby events. This implies that the local number density of double neutron star binaries 
or neutron star-black hole systems is significantly above previous estimates, if short GRBs 
come from such objects. In fact, it must be dominated by a so far undetected old population 
of double neutron stars or, alternatively, by observationally unconstrained NS+BH systems. 
A caveat of such conclusions and of comparisons to the observations is, of course, the still 
very small sample of well-localized short GRBs. Caution is therefore advisable until the 
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empirical foundation is more solid (see Bloom & Prochaska 2006). 

We restricted our analysis to the properties of the prompt emission of short GRBs, 
which may originate from neutrino-annihilation driven hydrodynamic jets, launched by the 
neutrino energy release of a hyperaccreting BH girded by a massive torus. Of course, one 
cannot exclude that magnetohydrodynamic effects might play a role that deserves more 
theoretical exploration. The discovery of X-ray flares following the short- duration, hard 
pulse of GRB 050724 after hundreds of seconds (Barthelmy et al. 2005) was interpreted 
as observational hint for the importance of such magnetic processes (Fan, Zhang, & Proga 
2005). If these late-time X-ray flares are indeed connected with an extended period of activity 
or a restart of the central engine long after the prompt emission is over (Burrows et al. 2005, 
Zhang et al. 2005, Liang et al. 2006), it seems not possible to explain these late outbursts by 
the accretion-rate dependent neutrino mechanism (Fan, Zhang, & Proga 2005). A prompt 
phase of neutrino (or neutrino plus magnetic) energy release that creates the short-hard 
GRB may thus be followed by a much more extended period of reduced source activity, in 
which X-ray flares could be powered by a magnetic mechanism (Fan et al. 2005). 

Admittedly, our investigation of the implications of the off-axis properties of short GRB 
jets from post-merger BH-torus systems was also based on a very limited set of numerical 
models (AJM05). This theory input of our analysis replaces the intrinsic peak luminosity 
function based on measured data used in other works (e.g., Nakar et al. 2005, Guetta & 
Piran 2005). The models were chosen to span approximately the possible diversity of the 
properties of ultrarelativistic outflows from such systems, but certainly represent it only 
incompletely. In particular, they were calculated with a single setup for the BH-torus system, 
in which the deposition rate of thermal energy and its spatial distribution (expressed in 
terms of the opening angle of the axial cone of the main energy deposition) was varied over 
a range guided by the results for neutrino-antineutrino annihilation in independent binary 
merger and post-merger evolution models. Neither was the torus evolved self-consistently 
in response to this energy release and to the jet formation, nor were magnetohydrodynamic 
effects included in the models. There is still a long way to go before the challenging problem 
of fully consistent modeling becomes manageable. This should include all relevant physics 
and track in full general relativity the history from the last orbits of the pre-merging binary, 
through the merging phase, to the possible BH formation and the subsequent secular post- 
merging evolution of the accreting relic black hole. Such simulations will ultimately have 
to be performed for a carefully selected sample of progenitor systems, which avoids the 
"arbitrariness" of the variations in our current set of models. Knowing that this goal is 
still far ahead and in view of the fact that the simulations available to us represent the 
present state-of-the-art of modeling collimated relativistic outflows from post-merger BH- 
torus systems, we tried to explore how far this information can help us bridging GRB engine 
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theory and observations. 
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Table 1: Observed short GRBs with determined redshifts. 



GRB a 


z 


Band 






Fluence 


E d . 

7,iso 


(L) e 


E P 






Gal. 


d h 


SFR 


Ref. 1 






keV 


ms 


ms 


10 — 8 erg cm -2 


10 50 erg 


10 50 ergs" 1 


keV 








kpc Mq yr 1 




050509b 


0.225 


15-150 


40 


33 


0.95 ± 0.25 


0.011 


0.8-1.5 


> 150 


1.5 ± 0.4 


2.3-4.3 


Ell 


40 


< 0.1 


1,5 


050709 


0.16 


30-400 


70 


60 


29 ±4 


0.167 


4.7 


83111 


0.7 ±0.2 


1.7 


SF 


3.8 


0.2 


2,6 


050724J 


0.257 


15-350 


250 
3000 


200 
2387 


11.3 ±8.7 
63 ± 10 


0.174 
0.97 


1.4-1.7 
0.7-0.8 


>350 


1.7 ±0.2 


1.6-1.9 


Ell 


2.6 


< 0.02 


3,7 


050813 


0.722 


15-350 


600 


350 


12.4 ±4.6 


1.6 


7.8-13 


>350 


1.2 ± 0.3 


1.7-2.8 


? 






4,8 


051221 


0.546 


20-2000 


1400 


906 


320±?'I 


24 


26 


400 ± 80 




1.0 


SF 


0.8 


1.5 


9 



a All the events where detected by Swift except GRB 050709, which was observed by HETE-2. 
b Observed duration Tqq. 

c Intrinsic duration (cosmologically corrected). 



d Isotropic equivalent gamma-ray energy released by the burst in the specified band (cosmologically corrected). 

e Average luminosity (computed as (L) = k E~ t i eo /T{). 

f Photon index of the average GRB spectrum (f(E) oc E~ v ). 

g Bolometric corrections to rest-frame energies of 20-2000 keV, following Bloom ct al. (2001). 
h Projected distance of GRB from host center. 

1 References: 1. Gehrels ct al. (2005); 2. Boer ct al. (2005); 3. Krimm ct al. (2005); 4. Sato et al. (2005); 

5. Bloom ct al. (2005); 6. Fox ct al. (2005); 7. Bcrgcr ct al. (2005); 8. Foley ct al. (2005); 9. Soderberg et al. (2006). 
J This GRB shows a first short pulse, followed by a softer and longer event. We report data for both events. 
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Fig. 1. — Properties of relativistic outflows from BH-torus systems as functions of viewing 
angle, determined by hydrodynamic simulations which show that collimated, ultrarelativistic 
mass ejection can be driven by sufficiently powerful deposition of thermal energy around the 
relic BH-torus systems of NS+NS/BH mergers, models B01-B08 differ in the rate of energy 
deposition per unit solid angle around the BH (AJM05). Left: Total (internal plus kinetic) 
isotropic equivalent energy of matter with terminal Lorentz factors > 100. This energy 
may be larger than the GRB energy by a factor depending on the efficiency of conversion of 
outflow energy to 7-rays. Right: Radially averaged Lorentz factor T at 0.5 seconds after the 
onset of energy deposition (AJM05). Note that the terminal value can be higher by a 
factor of 2-3. 
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Fig. 2. — Observationally relevant properties of relativistic outflows from BH-torus systems 
as functions of viewing angle, calculated using the results of the hydrodynamic models of 
Fig. I. Left: Total energy potentially radiated from ultrarelativistic, collimated ejecta with 
terminal Lorentz factors > 100, including the contributions from matter moving not along 
the line of sight (see text for details). The main difference compared to the left panel in Fig. 1 
is a smoothing of the jet wings and a slight widening at energies below about 10 49 erg. Middle: 
Relative observability of short GRBs from different viewing angles, estimated according to 
Eq. (2), using the -Ei so - values shown in the left panel. Right: Timescale of the observable 
GRB, estimated from the polar-angle dependent radial length A(#) of the ultrarelativistic 
outflow with Too > 100 in the observer frame, relative to the on-axis duration. 
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Fig. 3. — World-line diagrams for the evolution of shells at the forward and rear edges of 
the ultimately highly relativistic ejecta. While in the plot on the left the Lorentz factor is 
assumed to increase linearly and very quickly with distance from the source for both the 
leading and lagging shells, the plot on the right displays a situation that fits better the 
results of the hydrodynamic simulations of jet formation: The lagging shell accelerates more 
slowly than the leading one. The acceleration laws are indicated in the lower right corners 
of the two panels. Time t — is the moment when the shell at the rear edge of the ejecta 
is born, and A = ct cc is the radial extension of the outflow at this time, when t ce is the 
period of activity of the central engine. All lengths are normalized to A . In the figure on 
the left the radial thickness of the expanding shell nearly saturates at a value of about 2Ao, 
whereas on the right one can see a stretching of the fireball compared to its initial length by 
roughly a factor of 20, before a much slower growth indicates a nearly saturated state. The 
parameters for this plot were chosen to closely reproduce the results of one of the jet models 
of AJM05. 
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Fig. 4. — Expected normalized probability distributions of short GRBs as functions of red- 
shift z, fluence, isotropic-equivalent energy output in 7-rays, E 1} i so , and Lorentz factor T. 
The theoretical predictions are based on the results of our sample of jet models for E J: - lso 
and T shown in Fig. 2 and Fig. 1, respectively. In the upper panels we have adopted a 
lower fluence cutoff of 10 -8 ergcm -2 for the burst detectability, in the bottom panels we used 
1.6 x 10~ 7 ergcm~ 2 (note the corresponding change in the scales of the horizontal axes of the 
left two panels). The Lorentz factors T shown in the right panels are the values from the 
hydrodynamic jet models at 0.5 after jet formation. The terminal Lorentz factors Too will 
be significantly higher (see text). The lines correspond to three cases of different assumed 
comoving short GRB rate density as function of z (see text for details). The locations of 
the observed short GRBs listed in Table 1 are indicated. The very energetic GRB 051221 
is outside of the energy scale of the third panel and therefore it is omitted from all other 
panels of the figure as well. 



